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ABSTRACT 

We explore the effects of r-process nucleosynthesis on fall-back accretion in neu- 
tron star(NS)-NS and black hole-NS mergers, and the resulting implications for short- 
duration gamma-ray bursts (GRBs). Though dynamically important, the energy re- 
leased during the r-process is not yet taken into account in merger simulations. We 
use a nuclear reaction network to calculate the heating (due to /3-decays and nuclear 
fission) experienced by material on the marginally-bound orbits nominally responsible 
for late-time fall-back. Since matter with longer orbital periods i or b experiences lower 
densities, for longer periods of time, the total r-process heating rises rapidly with i or b, 
such that material with t or h ^> 1 seconds can become completely unbound. Thus, r- 
process heating fundamentally changes the canonical prediction of an uninterrupted 
power-law decline in the fall-back rate Mfb at late times. When the timescale for r- 
process to complete is > 1 second, the heating produces a complete cut-off in fall-back 
accretion after ~ 1 second; if robust, this would imply that fall-back accretion cannot 
explain the late-time X-ray flaring observed following some short GRBs. However, for 
a narrow, but physically plausible, range of parameters, fall-back accretion can resume 
after ~ 10 seconds, despite having been strongly suppressed for ~ 1 — 10 seconds after 
the merger. This suggests the intriguing possibility that the gap observed between the 
prompt and extended emission in short GRBs is a manifestation of r-process heating. 

Key words: nuclear reactions, nucleosynthesis, abundances-gamma rays: bursts 



1 INTRODUCTION 

One of the most important discoveries made with the Swift 
satellite is that short- and long-duration gamma-ray bursts 
(GRBs) originate from distinct stellar progenitors. While 
long duration GRBs track ongoing star formation (e.g. Fruc- 
ter et al. 2006) and result from the deaths of massive stars 
(e.g. Woosley & Bloom 2006), short-duration GRBs have 
now been localized to both early (Bloom et al. 2006; Berger 
et al. 2005) and late-type (Fox et al. 2005; Barthelmy et 
al. 2005) host galaxies, indicating a more evolved progenitor 
population (e.g. Prochaska et al. 2006). 

Although the origin of short GRBs remains unknown, 
the most popular and well-studied central engine model is 
the binary merger of two neutron stars (NS-NS; Eichler et 
al. 1989; Meszaros & Rees 1992; Narayan et al. 1992) or a 
NS and a black hole (NS-BH; Paczyhski 1991; Mochkovitch 
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et al. 1993) . This model is qualitatively consistent with both 
the host galaxy properties of short GRBs (Nakar et al. 2006) 
and the notable lack of a bright associated supernova in some 
well-studied cases (e.g. Hjorth et al. 2005). Depending on the 
detailed properties of the binary and the (unknown) supra- 
nuclear density equation of state, NS-NS/NS-BH mergers 
result in a central compact object (either a BH or hyper- 
massive NS) surrounded by a compact ~ 10~ 3 — O.IMq 
disk (e.g. Ruffert et al. 1996; Rosswog et al. 1999; Lee & 
Kluzniak 1999; Rosswog 2005; Faber et al. 2006; Shibata & 
Taniguchi 2006; see Lee & Ramirez-Ruiz 2007 and Faber 
et al. 2009 for recent reviews). The similarity between the 
estimated accretion timescale of this remnant torus and the 
typical duration of short GRBs (Tgo ~ 0.1 — 1 seconds) has 
also been used as evidence in favor of compact object merger 
models (Narayan et al. 1992). 

This clean picture has grown complex with the discov- 
ery that short GRBs are often followed by a "tail" of emis- 
sion (usually soft X-rays) starting ~ 10 seconds after the 
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GRB and lasting for ~ 30 — 100 seconds (Norris & Bon- 
nell 2006; Gelirels et al. 2006). While only ~ 1/4 of Swift 
short bursts show evidence for this extended emission, the 
observational limits are consistent with most bursts possess- 
ing an X-ray tail with a fluence comparable to that of the 
prompt GRB (Perley et al. 2009). Due to its rapid variabil- 
ity and other similarities to the prompt gamma-ray emission, 
the extended emission probably results from ongoing central 
engine activity. However, such a significant energy output 
on long timescales is difficult to explain in NS-NS/NS-BH 
merger models, most notably because the accretion disk is 
blown apart by a powerful outflow after only a few seconds 
of viscous evolution (Metzger, Piro, & Quataert 2008, 2009; 
Beloborodov 2008; Lee et al. 2009). 

One idea proposed for producing late-time activity is 
the "fall-back" of material ejected during the merger into 
highly eccentric (but gravitationally-bound) orbits (Ross- 
wog 2007; Lee & Ramirez- Ruiz 2007). If the ejected mass 
M has a distribution of energies dM/d\E\ oc |.E7| -a , mat- 
ter with E < produces a late-time accretion rate M 
onto the central object that decreases as a power-law in 
time: M oc (dM/dE)(d\E\/dt olh ) oc i -[( 5 + 2 «)/ 3 l (R ees 1988), 
where t OT h is the bound ejecta's orbital period or fall-back 
time. Indeed, a number of NS-NS/NS-BH merger calcula- 
tions find that the fall-back at late-times may be sufficient 
to explain the extended emission from short GRBs via con- 
tinued accretion (Faber et al. 2006; Rosswog 2007; Lee et 
al. 2009; see, however, Rossi & Begelman 2009). 

Fall-back accretion models rely on the assumption that 
matter ejected onto orbits with fall-back times t OI ^ 3> 1 
second remains bound. However, the gravitational binding 
energy of such long-period orbits is only 
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where t or b = 27r(a 3 /GM) 1//2 , a is the semi-major axis of 
the orbit, and M and m n are the central object mass and 
nucleon mass, respectively. 

The unbound ejecta from NS-NS/NS-BH mergers has 
long been considered a promising source for producing very 
heavy elements via rapid neutron capture (r-process) nu- 
cleosynthesis (Lattimer & Schramm 1974, 1976; Eichler et 
al. 1989; Freiburghaus et al. 1999). The total nuclear en- 
ergy available via the r-process (~ 1 — 3 MeV nucleon -1 ; 
see eq. [2]) greatly exceeds \E\ for orbits with t OT h J> 0.3 — 1 
second. As a result, r-process heating could have a crucial 
impact on the properties of late-time fall-back, an effect that 
has not yet been taken into account. 

In this paper, we examine the effects of r-process nucle- 
osynthesis on fall-back accretion in NS-NS/NS-BH mergers. 
In we describe the nucleosynthesis that occurs during 
the decompression from nuclear densities, and the extent to 
which NS-NS/NS-BH merger simulations properly capture 
the resulting energy release. This motivates fJ21 in which we 
present calculations of r-process heating along orbits that 
are nominally responsible for late-time fall-back. In 33] we 
discuss our results and their implications. 

2 DECOMPRESSION & ^-PROCESS HEATING 

Most of the material ejected when a NS is tidally disrupted 
originates from the NS's neutron-rich outer core, which has 
a typical electron fraction Y e ~ 0.1 set by 0— equilibrium 



under highly degenerate conditions (Pethick & Ravenhall 
1995; Haensel & Zdunik 1990a,b). Since the temperature 
remains fairly low as the ejecta expands (due to adiabatic 
losses), Y e probably remains low (~ 0.03 — 0.20) during the 
decompression from nuclear densities (Ruffert et al. 1997; 
Rosswog 2005). 

Schematically, the nucleosynthesis of decompressing 
neutron-rich matter can be divided into two stages: 

(1) Initial Decompression and Seed Formation (density 
P > Pdrip ~4x 10 11 g cm -3 ). During the earliest phases of 
decompression, very neutron-rich nuclei form, which rapidly 
emit neutrons as the material expands to lower density (Lat- 
timer et al. 1977; Meyer 1989). Heavy "seed" nuclei are 
then formed through (n, 7) reactions and, possibly, through 
charged-particle reactions in full nuclear statistical equi- 
librium (NSE). For example, Meyer (1989) finds seed nu- 
clei with average charges and masses Z ~ 40 — 70 and 
A ~ 90 — 110 (depending primarily on Y e and the expansion 
rate; cf. Goriely et al. 2004), while Freiburghaus et al. (1999) 
finds seeds with Z w 31 — 37 and A ~ 92 — 112. Since the seed 
nuclei and neutron mass fractions are given by X s — AY e /Z 
and X n = 1 — X 3 , respectively, the neutron mass fraction 
after initial decompression is large: X n ~ 0.3 — 0.9 for plau- 
sible ranges in the values of A, Z, and the electron fraction 
(Y e ~ 0.03 - 0.2; Ruffert et al. 1997; Rosswog 2005). 

(2) Rapid Neutron Capture (r-process) (p < pdrip)- Once 
the density decreases below neutron-drip, j3— decay channels 
begin opening in full, and a conventional r-process begins 
(see, e.g., Cowan, Thielemann, & Truran 1991 and Meyer 
1994 for reviews). In the r-process, very heavy nuclei (with 
peaks at A ~ 130 and 200) are formed when the seed nuclei 
rapidly capture the free neutrons remaining from stage 1. 
This establishes an (n,7) equilibrium, with f3— decays driving 
the nucleosynthetic "flow" to larger Z on longer timescales, 
with possible "fission cycling" between nuclei with A ~ 280 
and A ~ 130 — 140. 

If NSE is assumed, the nuclear energy released when 
seeds form (stage 1) can be captured in numerical simula- 
tions of NS-NS/NS-BH mergers by employing an appropri- 
ate equation of state (EOS). For instance, the NS-NS merger 
simulations of Rosswog et al. (1999) use a Lattimer-Swesty 
(1991) EOS, which accounts for the possible presence of pro- 
tons, neutrons, a— particles, and a single average "heavy" 
nucleus. As a result, they find that the ejected tidal tails "ex- 
plode" due to the energy released as seed nuclei form. Since 
the Shen EOS (Shen et al. 1998a,b) employed by Rosswog 
& Davies (2002; cf. Rosswog 2005) captures similar physics, 
the effects of seed nuclei formation are already taken into 
account in the fall-back estimates of Rosswog (2007). 

However, the subsequent r-process has the potential to 
generate a comparable or greater amount of energy. In par- 
ticular, once all of the synthesized nuclei decay back to stable 
isotopes, the total nuclear energy released and available to 
heat the ejecta is 



AE r ~ (1 - /„) 



x s 



X n A T 



(2) 



where A n = (m n — m p )c 2 = 1.293 MeV is the neutron-proton 
mass difference, and (-j) are appropriately-averaged 
binding energies for the seed and r-process nuclei, respec- 
tively. The factor /„ is the fraction of the nuclear energy 
lost to neutrinos and is ~ 0.5 (see £|3 - 2 p - Using typical val- 
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ues of (§) « 8 MeV nuc" 1 and (f ) « 8.7 MeV mic" 1 , 
we estimate that AE T ~ 1 — 3 MeV nuc - for X n in the 
range ~ 0.3 - 0.9. 

Comparing AE r with the binding energy of the fall- 
back material (eq. [I]), we conclude that if the r-process 
goes to completion, it will strongly affect the dynamics of 
orbits with t or b > 0.3 — 1 seconds. Unlike seed nuclei for- 
mation, the effects of r-process heating cannot be readily 
incorporated into merger simulations, in part because most 
of the energy is released on length and time scales exceeding 
that which can be presently simulated. More importantly, 
because the r-process is a non-equilibrium process involving 
a large number of exotic nuclei, its study requires a complex 
reaction network, which would be prohibitive to include in 
multi-dimensional simulations. In the next section we ex- 
plore the r-process heating of bound ejecta by performing 
nucleosynthesis calculations along a few representative La- 
grangian density trajectories. 

3 NUCLEOSYTHESIS CALCULATIONS 
3.1 Density Trajectories 

When the less massive NS is tidally disrupted during a NS- 
NS or NS-BH merger, a portion of the stellar material is 
ejected into one or two long tidal tail through the outer La- 
grange points (Lattimer & Schramm 1974). This material 
is imparted with a distribution of energies (or, equivalently 
for bound material, semi-major axes a). Initially, all of the 
ejecta (bound and unbound) is approximately spatially coin- 
cident (at an assumed pericenter distance r p « 10 7 cm) and 
shares a common density during decompression. Thus, dur- 
ing the early expansion we use the density trajectory p(t) 
corresponding to the unbound ejecta studied in r-process 
calculations by Freiburghaus et al. (1999) and taken from 
the NS-NS merger simulations of Rosswog et al. (1999). 

On later timescales, material with energy E < and 
fall-back time t or b oc | J E7 1 3 / 2 spatially decouples from the 
unbound ejecta, once their orbits approach apocenter. Mo- 
tivated by simulations (e.g. Rosswog 2007) and theoretical 
considerations (Rees 1988), we assume dM/d\E\ oc constant, 
corresponding to a fall-back rate of 

• (dM\ ( d\E\ \ 

~ 10- 2 M W dM/dE \( tmb Y% 

~ 10 M © S ^10-SMo/MeVnuc-iJ \U ) (3) 

where we normalize dM/dE so that ~ 10 _2 Mq returns to 
the central object on a timescale > 1 second (absent the 
effects of r-process heating), as is required in models that 
attribute late-time X-ray tails from short GRBs to fall-back 
(e.g. Rosswog 2007; Faber et al. 2006; Lee et al. 2009). 

The mass- flux along a series of trajectories with a given 
fall-back time torb can also be written as 

M = AQ.pvr 2 , (4) 

where r(t) and v(t) are the radius and velocity of the or- 
bit, and AQ(t) is the spread in solid angle of bound fluid 
elements. By equating equations <[3j and Q we obtain the 
late-time density trajectory. Our ignorance of the details of 
the merger (which depends on uncertainties such as the NS 
EOS) and the effects of nuclear energy input on the ejecta 
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Figure 1. r-process nucleosynthesis in bound ejecta from com- 
pact object mergers. (Upper Panel) Lagrangian density trajec- 
tories p(t) employed in our nucleosynthesis calculations, shown 
for ejecta with different initial orbital periods i or b = 0.1 s (short 
dashed lines), 0.3 s (dash — dot lines), 1 s (triple — dot — dash 
lines), 3 s (long dashed lines), and 10 s (solid lines) and for 
two values for the orbits' solid angle AQ = 4tt,4tc/10 (darker 
and lighter shaded lines, respectively). The trajectory of the 
unbound material (which all orbits share at early times) is shown 
with a dotted line. (Middle Panel) Temperature evolution for 
the trajectories shown in the upper panel. (Bottom Panel) Total 
r-process heating rate (due to j3— decays and nuclear fission) for 
the trajectories shown in the upper panel, assuming that 1/2 of 
the energy is lost to neutrinos (see !|3.2|l . 

trajectories are parameterized with AQ(t). For simplicity 
we assume that Afl is constant with time because our re- 
sults are relatively insensitive to this choice. We choose a 
relatively large value for Afl ~ 47r/10 — 47r, motivated by 
the large dispersion in the bound ejecta's orbital parameters 
expected to result from the explosive energy release during 
seed formation (Sj5J and the subsequent r-process heating. 
We only calculate heating until the orbits reach apocenter 
because we are primarily interested in the total energy re- 
lease, and r-process heating decreases rapidly once material 
re-compresses on its return to pericenter. The top panel of 
Figure [1] shows the density trajectories employed in our cal- 
culations for a variety of fall-back times. 

3.2 Network Calculations 

We use a dynamical r-process network calculation 
(Martinez-Pinedo 2008; Petermann et al. 2008) that includes 
neutron captures, photodissociations, f3— decays, and fission 
reactions. The latter includes contributions from neutron 
induced fission, /3 delayed fission, and spontaneous fission. 
All heating is self-consistently added to the entropy of the 
fluid following the procedure of Freiburghaus et al. (1999). 
The change of temperature is determined using the Timmes 



4 B. D. Metzger, A. Arcones, E. Quataert, G. Martinez- Pinedo 







1 








O 










' 



0.1 1.0 10.0 

U (s) 

Figure 2. Total r-process heating AE as a function of initial 
orbital period (or nominal fall-back time) t or b for the nucle- 
osynthesis calculations shown in Figure [T] shown for Af2 = 4n 
(asterisks) and 4-7r/10 (diamonds). The binding energy of the 
orbit \E\ (cq. [l]) is shown for comparison with a solid line (for 
an assumed M = 3Mq central object). 

equation of state (Timmes & Arnett 1999). Although our 
calculation does not explicitly account for the energy loss 
from j3— decays into escaping neutrinos, we take this into 
account by artificially decreasing the heating rate by a fac- 
tor 1/2. This is justified because most of the heating results 
from /3— decays and the energy released is shared approx- 
imately equally between electrons (which thermalize) and 
neutrinos. 

In addition to p(t), the initial temperature T, electron 
fraction Y e , and seed nuclei properties (A,Z) are specified 
for a given calculation. We assume an initial temperature 
T — 4 x 10 9 K, although the subsequent r-process heat- 
ing is not particularly sensitive to this choice (Meyer 1989; 
Freiburghaus et al. 1999). We also assume Y e — 0.1, Z ~ 36, 
A ~ 118 (see Varying the electron fraction and the 
properties of the seed nuclei will quantitatively affect the 
subsequent r-process heating (e.g. through the total avail- 
able energy; eq. [2]); the implications of this are discussed 
in JO] 

3.3 Results 

The middle and bottom panels of Figure [T] show the time 
evolution of the temperature T and heating rate E, respec- 
tively. The temperature initially decreases rapidly due to 
the adiabatic expansion. However, as the expansion rate 
decreases, nuclear energy generation reheats the ejecta to 
T ~ 5 x 10 8 K. At this point a difference develops between 
trajectories with short orbital periods (t rb < 1 s ) an d long 
orbital periods (i or b J> 1 s). For short orbital periods the 
density approaches an approximately constant value rela- 
tively quickly. At constant density all of the generated en- 
ergy is used to increase the temperature, which reaches a 
value T J> 10 9 K. At these high temperatures photodisso- 
ciation reactions become important and the r-process path 
moves closer to stability. This decreases the beta-decay rates 
and Q-values, and, consequently, the heating rate decreases. 
For the longest period orbits (i or b > 1 s), there is a balance 



between energy generation and adiabatic losses, which keeps 
the temperature almost constant. 

Also note that E remains relatively constant in the 
range ~ 2 — 5 MeV nuc -1 s _1 throughout most of the sam- 
pled trajectories for t < 1 second. This illustrates that the 

total r-process heating AE = J* * orb/ ' 2 Edt is dominated by 
relatively late times in the orbit, and that AE is approx- 
imately proportional to t OI ^. Note that at t ~ 0.7 second 
E sharply decreases once neutrons are exhausted and the 
r-process is effectively complete. 

The total r-process heating is shown explicitly in Figure 
[2] which plots AE for the trajectories from Figure [1] as a 
function of orbital period. The orbital binding energy \E\ oc 
*orb ( ec l- E) i s plotted for comparison with a solid line. 
Note that AE rises rapidly with i or b for both values of Afi, 
before saturating at AE r » 3 MeV (eq. [2]) for i or b > 2 
seconds; material with i or b > 1 second experiences sufficient 
heating to become unbound. 

4 DISCUSSION AND IMPLICATIONS 
4.1 Fall-Back Accretion 

Figure [2] illustrates that orbits with initial periods exceeding 
a fraction of a second experience sufficient r-process heat- 
ing to become unbound. However, this does not by itself 
guarantee a suppression in the fall-back rate Mo, at late 
times because the marginally-bound material ejected by the 
r-process (with initial orbital energy \Ei\ < AE) could in 
principle simply be replaced by material that was initially 
more tightly bound (\Ei\ J> AE). It is also important to un- 
derstand how the effects of r-process heating differ from that 
of seed nucleus formation; both release comparable amounts 
of the energy, yet seed formation is already incorporated in 
merger calculations and does not produce a sharp cut-off in 
the fall-back rate (Rosswog 2007). 

Although our r-process calculations assumed ballistic 
ejecta ( i)3.1|l . this approximation is no longer valid once sig- 
nificant energy is added to the orbit (i.e., if AE ~ \Ei\). 
Because the ejecta is optically-thick, most of the deposited 
thermal energy is transferred into kinetic energy via adia- 
batic expansion. This puts the ejecta on less-bound orbits 
which, due to their lower densities and longer periods, ex- 
perience even more r-process heating (Fig.[l|. This suggests 
that the r-process may lead to a run-away in which all ma- 
terial with an initial orbital period exceeding a threshold 
value will acquire sufficient energy to become unbound. 

To explore these issues quantitatively, we consider a toy 
model to calculate Mo, including the effects of r-process 
heating. We consider an ensemble of mass elements dM dis- 
tributed with initial energies Ei according to dM/d\Ei\ oc 
constant (as in eq. [3])- The final energy of each mass ele- 
ment at apocenter is given by 

/■*orb(B)/2 

E t = Ei+ I Edt, (5) 

Jo 

where E(t) is the heating rate along the orbit. As discussed 
in i )3.3l and shown in Figures [1] and [2] E is roughly constant 
in time along an orbit (for t < 1 s). Thus, we make the 

simplifying assumption that E = AE T /thea,t for t ^ thcat 
and E — for t > iheat, where AE r is the total available 
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Figure 3. Fall-back rate Mfb in NS-NS/NS-BH mergers including 
the effects of r-process heating, as calculated using the toy model 
described in i|4.1l Each model assumes a M = 3Mq central object 
and that the total available r-process energy is AE T = 3 MeV. 
Different line styles correspond to different heating timescales: 
thcat = 0-1 seconds (dotted line), 0.86 seconds (dashed line), and 
3 seconds (solid line). 



r-process energy (eq. [5]) and thcat is the timescale for re- 
process heating (which we leave as a free parameter). Our 
results presented below are relatively insensitive to the pre- 
cise functional form of E prior to thcat provided that the 
total heating is dominated by late times in the orbit. 

If the final energy at apocenter Ef is < for a given 
mass element, it remains bound despite r-process heating, 
with a new fall-back time which we approximate as t OI ^(Ef). 
If, on the other hand, Ef > the particle is unbound from 
the central object and does not contribute to late-time ac- 
cretion. Note that because the upper limit of integration in 
equation (0 increases with the orbital energy, this model 
allows for the run-away effect described above. 

Figure [3] shows our results for Mfb with r-process heat- 
ing, calculated for fixed AE T — 3 MeV and for several 
values of thcat- Because tightly-bound material with short 
orbital periods experiences little r-process heating, Mfb at 
early times is unaffected by the r-process and decreases at 
the canonical rate ex t~ 5//3 for all values of thcat- At late- 
times, however, there is a bifurcation in the behavior of 
Ma,: short heating times (thcat <C 0.9 s) lead to a relatively 
uninterrupted power law decline, while long heating times 
(thcat > 0.9 s) produce a sharp cut-off in the fall-back rate. 

The origin of this bifurcation can be understood by not- 
ing the existence of a critical orbit distinguishing material 
that necessarily experiences the full heating available AE r 
from those that may not. When r-process heating is active, 
the energy of an orbit evolves as E = — \Ei \ + Et. As time 
increases, the magnitude of the orbital energy decreases, i.e., 
\E\ goes down, and thus the orbital period increases. There 
is a critical orbit for which the r-process heating leads to the 
orbital period increasing so rapidly with time that the orbit 
can never actually reach apocenter (t = t or b/2) as long as 
r-process heating is active. These orbits necessarily receive 
the full r-process heating AE r . This critical orbital energy 
E c can be determined by setting dt/d\E\ — l/2(dt OI b/d\E\) 
using equation (fTJ and E = — \Ei\ + Et. This implies 



The initial orbital energy for which E c is just reached 
at apocenter is given by E(t — t rb/2) = E\ + Et OI b/2 — E c , 
which implies \E- ljC \ — 5\E c \/3. Note that equation ((6J can 
be estimated on dimensional grounds by solving for Ei such 
that \Ei\ » {AE r /t hcat )t olh (Ei)/2. 

There is a bifurcation in the behavior for \E{\ < \Ei ;C \ 
and \Ei\ > |£?i lC |. Orbits with < |i?i, c | necessarily receive 
the full heating and so have Ei — — \Ei \ + AE T , while those 
with \Ei\ > |-Ei, c | may or may not (see below). The orbital 
period (or fall-back time) corresponding to |i? c | is given by 




from which it follows that 



i*22l ~ 1 7 r I hcat | (f(\ 

W ' V 3M 0/ V 3 MeV J V Is ) ■ K > 

If theat <C torb,c then any material with t or b > t or b.c has 
already experienced the full heating AE r earlier in its orbit. 
In this case, there is a slight decrease in Mfb around the time 
at which t or b ~ theat , but there is no significant interruption 
in the fall-back rate; this corresponds to thcat = 0.1 s in 
Figure [3] In particular, Mfb still decreases as a power-low oc 
t -5 / 3 at late times because adding a constant energy to each 
dM simply renormalizes the energy scale when dM/dE is 
flat. This explains why seed nucleus formation, which occurs 
on roughly the initial expansion timescale ~ milliseconds 
<C torb.c, has little effect on the rate that Mfb decreases at 
late times. 

On the other hand, if thcat 2> t or b,c then an absolute 
cut-off in dM/d\E{\ (and hence Mfb) occurs for Ei :> E c , 
corresponding to times t ^> torb.c- This case is well-illustrated 
by the thcat = 3 s model in Figure [3] 

In intermediate cases, when thcat ~ torb.c (i.e., |-Ei, c | ~ 
AE r ), there is still cut-off in the accretion, but material with 
\Ei\ > \Ei, c \ may remain marginally- bound despite the extra 
energy it receives, thus leading to a temporal gap in Mfb(t). 
This is illustrated by the thcat = 0.86 s model in Figure [3] 
which shows a long delay between the cut-off in accretion at 
t ~ 0.4 s and its resumption at t ~ 10 s. We find, however, 
that the gap only exists for a fairly narrow range of param- 
eters, and, when present, its width At gap is exponentially 
sensitive to t hcat /t OThtC ; e.g., increasing thcat/t or b,c from 1.3 
to 1.7 increases At gap from ~ 1 to ~ 100 s. 

From our calculations in ^3]we find that AE? ~ 3 MeV 
and thcat ~ 0.7 s, which corresponds to thcat ~ l-5t orb ,c- 
This is at the boundary between the "absolute cut-off" and 
"intermediate" regimes described above and shown in Fig- 
ure We discuss the implications of this result in the next 
section. 

4.2 Implications for the Origin of Short GRBs 
with Extended Emission 

Our primary conclusion is that fall-back accretion following 
NS-NS/NS-BH mergers is suppressed on timescales exceed- 
ing ~ 0.3 — 1 seconds due to r-process heating. This result 
has important implications for the origin of extended X-ray 
emission observed ~ 10 — 100 seconds following some short 
GRBs. As discussed in the Introduction, standard merger 
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Figure 4. Fall-back rate Mft, calculated using the model from 
i]4.1l taking E(t) directly from our r-process calculations in Figure 
[T]for the t or b = 10 s orbit. 



models have difficulty explaining activity on such a late 
timescale, which has lead to the suggestion that late-time 
flaring is powered by fall-back (Faber et al. 2006; Rosswog 
2007; Lee et al. 2009). If long heating timescales and/or high 
values of AE r obtain (such that ihcat > 2t or b,c), our results 
strongly disfavor this explanation because the cut-off in Mfb 
is absolute: very little material returns to the central object 
at late times (see the solid line in Fig. [3} . 

This conclusion changes, however, if AE r is lower 
and/or thoat is shorter, such that theat < 2i or b, c - In this case, 
the r-process produces a temporal gap in the fallback rate 
instead of an absolute cut-off (Fig. [3]). Intriguingly, the ex- 
tended emission following short GRBs shows a lull of ~ 3—10 
seconds between the end of the GRB and the beginning of 
the extended emission (e.g. Norris & Bonnell 2006; Gehrels 
et al. 2006; Perley et al. 2009). Attributing this delay to r- 
process heating appears, however, to require fine-tuned pa- 
rameters: theatAorb.c must be between ~ 1.5 — 1.7 in order 
for the temporal gap to have a duration of 3 — 30 seconds. 
Nevertheless, typical parameters for r-process heating are 
not far from the critical condition i ne at ~ iorb.c- 

To illustrate this explicitly, Figure [3] shows Mfb(f) cal- 
culated using the simple model described in i|4.1l but using 
E directly from our r-process calculations in Figure Q] for 
iorb = 10 seconds. For these particular parameters, the ex- 
act E(t) suppresses fall-back accretion from ~ 0.4 — 4 sec- 
onds, but fall-back then resumes at late times. Although 
the properties of the gap are quite sensitive to the param- 
eters of the r-process heating, the presence of a gap may 
in fact be rather robust, because the parameters on which 
thcat and t orhtC depend (AE T , M,t hBElt ) may not vary sub- 
stantially from event to event. For instance, the outer por- 
tions of a NS ejected during a NS-NS/NS-BH merger are 
probably highly neutron-rich and likely remain so during 
the expansion (Ruffert et al. 1997). For low Y e the total 
energy released by the r-process asymptotes to the value 
AE r ~ 0.5[(B/A) T - A n ] ps 3.3 MeV (provided that Y e is 
not sufficiently low that the r-process freezes-out before all 
neutrons are captured). Furthermore, for a merger to pro- 
duce a significant accretion disk (as is required to produce 
a GRB), the mass of the central black hole or NS proba- 



bly must lie in the relatively narrow range M ~ 3 — IOMq 
(e.g. Rantsiou et al. 2009). Finally, a heating timescale of 
~ 1 second also appears to be robust, relatively indepen- 
dent of uncertainties such as Afl (see Fig. [2}; £ nC at ~ 1 sec- 
ond was also found in an independent r-process calculation 
employing somewhat different assumptions and physics (see 
Goriely et al. 2004; their Fig. 3). Further work is thus clearly 
required to assess the possibility that the gap between the 
prompt and extended emission in short GRBs is a manifesta- 
tion of r-process heating. In particular, given the sensitivity 
of the late-time fall-back rate to the r-process heating, it is 
important to consider a range of initial Y e and seed nuclei 
in the r-process calculations, and to improve on the simple 
model in i]4.1l by carrying out hydrodynamic calculations 
of fall-back using simplified models of the r-process heating 
rate. 

If future work indicates that r-process heating produces 
a robust sharp cut-off in Mfb (with no late-time resumption) , 
our results may instead suggest that the central engine in 
some short GRBs is a long-lived NS rather than a BH, be- 
cause the former could remain active even in the absence of 
surrounding matter. For instance, the massive central object 
that forms in a NS-NS merger could be temporarily sup- 
ported by differential rotation (e.g., Baumgarte et al. 2000; 
Duez et al. 2004, 2006) or remain stable indefinitely if it 
loses sufficient mass via a centrifugally-driven outflow (e.g. 
Thompson et al. 2004; Dessart et al. 2008). Another possi- 
bility for producing a stable NS is via the accretion-induced 
collapse (AIC) of a white dwarf. In either AIC or NS-NS 
mergers, if the rapidly-rotating NS is strongly magnetized, 
its electromagnetic spin-down could plausibly power the ob- 
served extended emission (Usov 1992; Metzger, Quataert, 
& Thompson 2008). If AIC occurs following a double white 
dwarf merger, late-time emission could also be powered by 
the accretion of material left over after the merger, which has 
an accretion timescale ~ 100 seconds (Metzger, Quataert, & 
Thompson 2008). 
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